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Abstract. We successfully utilise the SPARC code to model fast-Alfven mode conversion in 
the region ca 3> cs via 3-D MHD numerical simulations of helioseismic waves within constant 
inclined magnetic field configurations. This was achieved only after empirically modifying the 
background density and gravitational stratifications in the upper layers of our computational 
box, as opposed to imposing a traditional Lorentz Force limiter, to ensure a manageable 
timestep. We found that the latter approach inhibits the fast-Alfven mode conversion process 
by severely damping the magnetic flux above the surface. 



1. Introduction 

A series of recent studies (e.g., [Il[2l[3l[Il[5j|6l[71[8]) have shown that an important feature not 
widely accounted for in sunspot seismology is fast-Alfven wave mode conversion - a process which 
occurs at, and beyond, the fast wave reflection height (where ca ~ u / k^; ca denotes the Alfven 
wave speed, oj the wave frequency and the horizontal wavenumber) - being spread over many 
scale heights for wavenumbers typical of local helioseismology. This process is most efficient 
for 9 (field inclination from vertical) between 30° — 40°, and <j> (angle between the magnetic 
field and wave propagation planes) between 60° — 80°, and appears to have the potential to 
modify the seismic wave-path through the solar atmosphere, thereby affecting the wave travel 
times that are the basis of our inferences about the subsurface (see [9j [TQl HH [12] for recent 
reviews). Motivated by these studies, our aim is to use the Seismic Propagation through Active 
Regions and Convection (SPARC) code, a 3-D magnetohydrodynamic (MHD) wave-propagation 
code developed by |13| for computational heliosiesmology, to numerically simulate this process, 
and investigate the implications of fast-Alfven wave mode conversion on the seismology of the 
photosphere. In this paper we describe our attempts at forward modelling this process with 
SPARC using a quiet-Sun background model permeated by homogenous inclined magnetic fields. 

2. Numerical Setup 

The SPARC code solves the 3-D linearized Euler and induction equations of magnetofluid motion 
in Cartesian geometries to investigate wave interactions with local perturbations (e.g., sound 
speed, pressure, density, flows, magnetic field etc.). Over the past few years, a number of various 
solar phenomena have been studied using this code (e.g., [21 [151 EE EH ttH] ) • The computational 
box we employ for our simulations using SPARC spans 186.6 Mm in the horizontal direction 
(128 evenly-spaced grid points in the horizontal directions x and y; Ax = Ay = 1.46 Mm/pixel), 



and from 2.5 Mm above the surface (z = 0; where z denotes height in Mm) to 25 Mm below in 
the vertical direction (300 non-uniformly spaced grid points in z; Az varies from several hundred 
kilometres at depth, to tens of kilometres in the near-surface layers). The vertical boundaries 
of the box are absorbent, with perfectly matched layer (PML) boundary layers spanning the 
top 10 and bottom 7 grid points in z, while periodic boundary conditions are imposed on the 
horizontal sides. In a similar manner to El [7], we use a monochromatic [y = uj/2tt = 5 mHz) 
plane-parallel wave driver, which only excites waves propagating in the (x, z) planes. We do 
this by imposing a perturbation of the form: 

sin(wt) e -(x-x°) 2 /z% & -{z-z a ?/^l (1) 

in a few grid points near zq = — 5 Mm and centred at around xq = Mm, in pressure, density 
and velocity (x and z components only). 

3. Model Atmosphere 

The background model used is a convectively stabilised solar model (CSM_B) from |19j . On top 
of this background model we employ a constant, inclined magnetic field configuration using the 
prescription from pQ: 

Bo = Bo (sin cos 0, sin 9 sinc/>, cos <p), (2) 

We choose Bq = 1500 G and a number of different angle orientations in the range < 6 < 90° 
and < (j)< 90°. 

In 3-D MHD simulations, the timestep (At ~ Az/ca) is often highly constrained by 
the Courant-Friedrichs-Lewy (CFL) condition, due to the exponentially increasing value of 
ca = Bo/(/Uo/o) 1//2 (where p denotes density and po the magnetic permeability) above the surface. 
This causes the wavelengths of both the fast and Alfven waves to become quite large, resulting 
in an extremely stiff numerical problem. The most common way of dealing with this issue, in 
computational helioseismology, has been to simply apply a ca "limiter" to moderate the action 
of the Lorentz Force when the ratio ca/cs (where cs denotes sound speed) becomes exceedingly 
large. Some choices for these limiters have been discussed before (e.g., [201 ED 113 1221 
Generally, they tend to prefix the Lorentz-Force terms in the momentum equations, with ca 
typically being capped at ~ 20 — 60 km/s. The physical implications of artificially limiting 
the Lorentz Force in such a manner (particularly on the seismology) have not been explored. 
Recently though [24J has shown that significant internal reflection of Alfven waves can occur if 
the ca profile is not adequately treated above the surface. 

Our method of ensuring a reasonable At for our calculations involves empirically modifying 
the background density (po(z)) in the upper layers (~ 0.5 < z < 2.5 Mm) of CSMJB, in 
conjunction with a commensurate modification of the gravitational acceleration (go(z)) profile 
over the same z range (in order to offset the change in the density scale height), to obtain 
a maximum global fast (cp) and ca of ~ 80 km/s. All other background variables remain 
unaltered. The modified profiles of po(z) and go(z) are shown in Figure [lj The value of 80 km/s 
was chosen because it provides us with a reasonable time step (At = 0.5 s), and is safely higher 
than the largest horizontal phase-speeds (iv/kh) typically sampled in sunspot seismology (e.g., 
[25]). This is important since ui/kh ~ ca also denotes the location of the fast mode reflection 
height in the solar atmosphere [6]. Figure [2] a) shows the resulting ca and cp profiles in our 
model as a function of height. For comparison purposes, the ca and cp profiles which would 
result from imposing a Lorentz Force limiter, instead of modifying po(z) and go{z) above the 
surface, are shown in Figure [2] b). 

The obvious downside of empirically modifying po(z) and go(z) in order to satisfy the 
CFL condition is that background model will no longer be as 'solar-like' (i.e., in terms of 
eigenfrequencies, eigenfunctions and power spectrum) as CSM_B. The larger po(z) profile which 
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Figure 1. Plots of the original CSM_B ([19j: dot-dashed line) and modified (solid line) 
background density (left) and gravitational acceleration (right) profiles as a function of height. 
The solid vertical line represents the ca ~ cs height. 



now results above the surface also modifies acoustic cut-off frequency (u ac = cs/2H p ; where H p 
denotes the density scale height), which is reduced from v = 5.2 mHz to 3.6 mHz. We also find 
that the modified atmosphere produces large-amplitude convective (g-) modes at v « 1.5 — 1.7 
mHz (it is worth noting though that this is a frequency range which is typically associated with 
supergrannulation noise and is generally filtered out /ignored in sunspot seismology). 
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Figure 2. a): Plots of various wave speeds as a function of height resulting from the modified 
CSM_B background model, b): Same as a) but for an unmodified CSM_B background in 
conjunction with a Lorentz Force limiter (i.e, the Lorentz Force is multiplied by a factor of 
200c|/(200c| + c A ), leading to ca and cf being capped at « 80 km/s). The solid vertical line 
represents the ca ~ cs height. 



However, more importantly for our concerns, this method ensures that we satisfy our CFL 
condition without any direct modification of the Lorentz Force via the introduction of an 



artificial term in Maxwell's equations, which as we shall show in the proceeding section, results 
in unphysical damping of the magnetic flux above the surface and inhibits the fast- Alfven mode 
conversion process. 



4. Results 

Following [5j (H |7|, we use velocity projections onto three orthogonal directions (e; ons , which 
selects the longitudinal component of wave propagation, i.e, the slow mode; e trans , which selects 
the transversal component of wave propagation, i.e., the fast mode; e perp , which selects the 
perpendicular component of wave motion, i.e., the Alfven mode; see equations 1-3 in [5] for 
definitions) to separate the Alfven mode from the fast and slow magneto-acoustic modes in 
the region ca S> C5. We also calculate the temporally averaged acoustic (F ac = (piVi); where 
p represents pressure and v represents the 3-D velocity respectively, subscript "1" represents 
perturbations), magnetic (F mag = (Bi x (vi x B )//io)) and total (F to t = F ac + F mag ) energy 
fluxes in order to measure the efficiency of conversion to Alfven waves around the ca ~ cs 
equipartition height. 




Figure 3. Figures above represent the velocity projections derived from simulations using 
constant inclined magnetic field configurations (left column: 6 = 0°,(p = 80°; middle column: 
9 = 30°, 4> = 80°; right column: 9 = 60°,</> = 80°) as a function of height (Mm) and time 
(minutes). V perp and Vtrans amplitudes have been scaled by a factor of (pqca) 1 ^ 2 , while Vi ong 
has been scaled by a factor of (pocs) 1 ^ 2 - The grayscale is the same in all panels. The horizontal 
solid and dotted lines represent the ca ~ cs and fast-mode reflection heights respectively. 



The results of the projected velocities derived from simulations where 4> is fixed at 80° and 
9 varies from 0° to 30° to 60° are shown in Figure [3j The inclination of the ridges in these 
projections indicates the wave propagation speed: the more inclined the ridges, the lower the 
propagation speed and vice versa. The presence of the (primarily acoustic) slow mode (Vi ong ) is 
clearly visible above the ca ~ cs level (solid line) in all three cases (with the reduced uj ac resulting 
in a significant amount of acoustic waves propagating along the field both above and below the 
equipartition height), while the rapidly propagating Alfven mode (y perp ) only appears when the 
magnetic field is sufficiently inclined and oriented out of the plane, i.e, 9 = 30°, (ft = 80° and 
9 = 60°, <p = 80°, with the latter configuration appearing to be the more efficient in producing 
Alfven waves. For these two cases, a faint presence of the magnetically dominated fast mode 



(Vtrans) above the reflection height (dotted line) can still be made out. This is a result of 
our plane-parallel driver, which excites all wavenumbers, leading to a proportion of fast modes 
being transmitted through to the PML, rather than being reflected at ca ~ cj/kh. In the 
9 = 60°, <j) = 80° projection, there appears to be a wavefront with an opposite inclination of 
ridges close to the top boundary of the domain. This could either be an artefact of the colour 
scheme/scaling, or a numerical artefact (i.e, reflection) from the upper boundary condition 
of the simulations. While we do not observe any reflection from the upper boundary in the 
corresponding acoustic flux (see Figure [4|, given that the wavefornt appears to arrive at the 
upper boundary prior to the arrival of the slow waves, it could also be possible that the signature 
is a yet unmodelled product of the fast-Alfven conversion process. This is something which we 
hope investigate in a future work. 

Figure [4] shows the vertical component of the averaged fluxes as a function of height. We 
observe that the flux variations are strongest near the conversion layer (solid vertical line), with 
the magnetic flux exceeding the acoustic flux when 9 = 60° , <fi = 80° for z > 1 Mm. These results 
are in good agreement with previous numerical simulations of fast-Alfven mode conversion using 
homogenous inclined magnetic fields (e.g., see Figure 3 from [5]). 
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Figure 4. Time-averaged magnetic (bold sold line), acoustic (dashed line) and total (dotted 
line) fluxes (in non-dimensional units) calculated as a function of height. The solid vertical line 
represents the ca ~ cs height. 



For comparison purposes, we also conducted simulations where, instead of modifying po(z) 
and go(z) above the surface, we employ a Lorentz Force limiter (with ca and cp capped at 80 
km/s, i.e, as shown in in Figure [2] b)). The resulting averaged fluxes are shown in Figure [5j 
While differences in the magnitude and height variations of the acoustic fluxes between these 
results, and those contained in Figure |4j can be explained by the change in u ac , the differences 
in the magnetic fluxes, particularly when considering the 9 = 60°,</> = 80° cases, are almost 
entirely due to the Lorentz Force limiter. With the Lorentz Force limiter in place, the magnetic 
flux above z > 1 Mm appears to just be able to creep above the acoustic flux for a couple 
of hundred kilometres, before being completely damped prior to reaching PML. We observed 
this phenomenon regardless of the value of the ca cap that was used with the Lorentz Force 
limiter. As expected, the resulting velocity projections for these cases (figures not included) also 
confirmed the absence of any significant Alfven modes above ca ~ cs- 



5. Summary 

Understanding the physics of propagating waves within regions of strong magnetic fields, of 
which fast-Alfven wave mode conversion has recently been shown to be a critical component, is 
essential for helioseismic studies of sunspots and active regions. We used the 3-D linear MHD 
solver SPARC to simulate this process in a convectively stabilised solar model (CSM_B with 
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Figure 5. Time-averaged magnetic (bold sold line), acoustic (dashed line) and total (dotted 
line) fluxes (in non-dimensional units) calculated as a function of height for simulations where a 
Lorentz Force limiter is used with ca and cp capped as shown in Figure [2] b). The solid vertical 
line represents the ca ~ cs height. 



empirically modified po(z) and go (z) profiles in the upper layers to ensure a reasonable At) 
permeated by homogenous inclined magnetic fields. We found that employing a traditional 
Lorentz Force limiter to artificially cap ca above the surface tends to inhibit the fast-Alfven 
mode conversion process by significantly damping the magnetic flux above the surface. 

The next steps in our forward modelling process will include the introduction of random 
stochastic sources and more realistic (i.e., sunspot-like) background atmospheres, in order to 
simulate artificial helioseismology data sets. With the aid of local helioseismic diagnostic tools, 
such as time-distance helioseismology and helioseismic holography, we will then be able to 
attempt to quantify the effects of fast-Alfven mode conversion on the wave travel times. 
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